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ABSTRACT 


Observations of the 3^3 - 230 transition of water vapour in the 
direction of Ori MCI in 1980 February show a 50% flux increase and an 
apparent additional red shift of approximately 2 km s' 1 relative to the 
line observed in 1977 December. We have adapted the mm-wave telescope 
calibration formalism to suit optical telescopes, and included an explicit 
relation to aperture efficiency to facilitate the use of extended source 
calibration data with spatially unresolved sources. From a detailed 
examination of the amplitude and frequency calibration, it appears unlikely 
that the effect is due to systematic error. The increase is attributed to 
the appearance of a new component at a velocity of 12 km s~* with respect 
to the local standard of rest. The new component also has broad wings. 
Increased emission from a region in the high-velocity core of Ori MCI can 
be due either to additional far-IR radiation to pump the 183 GHz transition 
or to a change in the physical conditions in the gas. We have carried out 
statistical equilibrium calculations using the large-velocity-gradient 
formalism to develop a model for the emission and to identify the physical 
conditions to which the excitation of the 3 j3 - 2 q 2 transition is sensitive. 
The calculations support a model in which the gas in the region of enhanced 
emission is hotter than the dust. The temporal coincidence between the 
183 GHz increase and the 22 GHz water maser outburst suggests a common, 
impulsive cause, which has heated the gas in a part of the HV source, 
enhancing the emission in both transitions. 

Subject headings : interstellar: molecules — masers — stars: 

formation — nebulae: Orion — radio sources: lines 
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I. INTRODUCTION 


Among the features which make the core of Ori MCI a unique source of 
molecular line radiation are strong emissions at 183 GHz due to the 3j^ - 
220 transition of water vapour (Waters et_ aJ_. 1980a) and at 380 GHz due to 
its 4^4 - 321 transition (Phillips, Kwan, and Huggins 1980). Both lines 
have a shape similar to many molecular lines, consisting of a narrow 
central "spike" surrounded by broad Doppler wings which are referred to as 
a "plateau" (Zuckerman and Palmer 1973). The wings of the 183 GHz line were 
much stronger than had been observed until that time in other molecules. 
Assuming that the water vapour wings came from a region at least as small 
as that which gives rise to the wings in CO, it was found that a mild 
population inversion was required in the 3^3 - 2 q 2 transition to explain 
the observed intensity (Kuiper, Zuckerman, and Rodriguez Kuiper 1981). 

In 1980 February, brief observations at 183 GHz were made of Ori MCI 
as an instrumental check while conducting a search for 183 GHz emission 
from late type stars (Kuiper et al . 1984). Relative to the previous 
observations of 1976 and 1977, the peak antenna temperature was nearly 
twice as strong. We report these results and consider the extent to which 
the previous interpretations need to be modified to account for this 
change. 

Water maser emission at 22 GHz is a recognized signpost of star 
formation. Yet, an understanding of the conditions which give rise to this 
emission continues to be elusive (e.g., Elitzur 1982). We consider the 
implications of our observation for the 22 GHz maser outburst which was 
observed during the same period. 
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II. OBSERVATIONS 


The observations were made with Schottky diode mixer receivers 
mounted at the bent-Cassegrain focus of the 91 cm telescope of the G. P. 
Kuiper Airborne Observatory. The instrumentation for the 1977 observations 
has been described by Waters et al. (1980a, b). A cooled receiver was used 
for the observations of Orion in 1980 February and is described in 
subsection a below. 

During the 1977 observations, the calibration measurement consisted of 
an observation of the moon. The moon was not available in the 1980 February 
flight series, but was observed during a subsequent flight series in 1980 
August (when Orion could not be observed). However, since the telescope 
did not change at all between these flight series, and the receiver only 
minimally (see below), we feel that the 1980 August calibrations can be 
applied to the 1980 February data. 

In addition, in 1980 August, we made measurements inside the telescope 
structure and of the cold sky which, when combined with calculations based 
on our measured horn power pattern, identified 85% of the losses in an 
absolute sense. Since the telescope structure did not change from 1976 to 
1980, we feel that the 1977 and 1980 observations of Orion can be compared 
to an uncertainty of ~ 10%. The lunar observations confirm this for those 
efficiency factors which are not related to the beamwidth. The remaining 
factor for unresolved sources depends on aperture illumination which cannot 
have departed from the designed pattern sufficiently to account for an 80% 
change in antenna temperature. 

Readers who do not wish to study all the details are invited to skip 
to subsection c. Results. 
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a) Instrumentation 


In 1980 February, calibration was accomplished by intercepting the beam to 
the telescope with a flat plate and diverting it downward to a rotating 
semi-circular mirror which directed it alternately to an ambient absorber 
and an absorber immersed in liquid nitrogen. To verify the operation of 
the spectrometers by measuring the atmospheric 183 GHz emission, we used a 
different configuration, which we describe below, to make our calibration 
measurements. However, except for the difference in mounting the receiver 
to the telescope, the systems were identical. 

A simplified diagram of the receiver, the KAO telescope, and optical 
path (not to scale) are also shown in Figure 1. The detector was a Schottky 
diode, mounted in a reduced height waveguide and cooled to approximately 
15 K by a closed-cycle helium refrigerator. Power was coupled from the 
waveguide mode to the Gaussian beam mode by a small corrugated circular 
horn and a Teflon lens. The local oscillator .power was injected through 
a quasi-opt ical , folded Fabry-Perot diplexer. The intermediate frequency 
was amplified by a gallium arsenide field-effect transistor amplifier, 
cooled to approximately 80 K. The fixed oscillators used in the frequency 
conversions were counted at JPL before or after each flight series and 
found not to differ from their nominal frequencies by more than 2 kHz. The 
data system consists of a 36-channel filter bank with resolution at the 
center frequency of 1 MHz tapering to 8 MHz at the edges and a 256-channel 
digital Fourier transform spectrometer with 40 KHz resolution over a 10 MHz 
bandwidth. 

The antenna temperature scale was calibrated relative to the ambient 
and liquid-nitrogen cooled targets in the receiver. The temperatures of 
the load were measured by means of attached thermistors. The receiver's 
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internal calibration was checked in 1980 August by placing ambient and 
liquid-nitrogen soaked Eccosorb foam targets at the receiver aperture. 

We concluded that the receiver's internal calibration was accurate to 1% 
in slope, and 0.3% in zero offset, relative to the thermistor temperature 
readings. The double sideband system noise temperature was typically 
900 K for the operational system. 

The internal optical components of the receiver were aligned by shining 
a laser down the optical axis, replacing the opaque Teflon lenses by small 
irises. The receiver was aligned with the KAO optics by focusing the image 
of Polaris on the center of the pressure window. Because of the long focal 
length and long wavelength, the Schottky diode detector is required to be 
only within 2 mm of the optic axis and several cm from the focus along the 
optic axis. We believe that the pointing was accurate to ~ 10 arcsec. 
Pointing was verified during both the 1977 and 1980 flights by making a 
five-point map of Orion, which showed that our. pointing error was less than 
1/2 arcmin. 

To remove sky and telescope background, beam switching between source 
and adjacent sky was used, subtracting sky from source in real time. The 
relative positions of the two beams were reversed in alternate integrations 
which were then subtracted to remove baseline effects due to asymmetry 
between the two beams. Any residual baseline effects were due to asymmetry 
in the system and variations on the timescale of several minutes. 

b) Telescope Calibration 

We consider how to apply the overall efficiencies determined from 
observations of the moon in 1976 October and 1980 August to observations 
of unresolved sources. We consider both the losses in the telescope and 
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through the atmosphere. 

Following the notation of Kutner and Ulich (1981; hereafter referred to 
as KU), the change in antenna temperature due to the presence of a source 
may, for single sideband observations, be written as 

— r a (v)A 


AT A “ n l ^fss "c T R e 3 


( 1 ) 


The p's are efficiency terms, where 0 < n < 1; Pi accounts tor losses in the 
telescope that terminate in the ambient temperature surroundings; Pf ss is 
the efficiency due to spillover past the secondary (we assume that there is 
no diffuse scatter from optical reflecting surfaces at our wavelength); p c 
is the efficiency with which the diffraction beam couples to the source; 

T a (v) is the zenith atmospheric opacity at frequency v; A is the airmass; 
and Tjj is the radiation temperature of the source. The telescope efficiency, 
rij, can for convenience be divided into three terms 

hi = h r Hint h mo (^) 

where p r is the efficiency due to ohmic losses in the telescope's reflecting 
surfaces; Pint accounts for losses between the receiver and a point adjacent 
to the tertiary flat blocking the secondary; p mo is the efficiency of the 
main optics accounting for blockage of the primary by the secondary and the 
spider. We ignored blockage by the spider. (If the integral over the 
"forward hemisphere" in KU is understood to include only the part of the 
power pattern which sees the sky, then q^ nt P m0 is P rss in the notation of 
KU.) 

When the source is small relative to the telescope beam, Equation 1 
may be rewritten (see Eq. 3 of KU) as 


AT A ’ = AT A e 


T a ( v)A 


T R // 

// p n (a)dn “ s 

4ir 
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hr A 2 S 

= = n A AS/2k 

//, P n (W)dO 2k 
4ir 

where S is the flux density and the aperture efficiency, p A , has its usual 
definition: 

A 2 /A 

p A = n r — n r A 2 /A% 

// P n (fl)dQ 

4 11 __(3) 

Because n A changed between observations, we compare a modified antenna 
temperature 

A T A 

T A /n A e 3 T R // B n (f)d» 

a a x2 K '' ns n 

-x A 

= AS e a /2k (4) 

(We note in passing that n A = p g A 2 /A%, where n s = Pf ss Pj , the extended 
source efficiency of KU, is referred to in early literature as the beam 
efficiency. ) 

During the flight series in 1980 August, we placed ambient temperature 
and liquid nitrogen-cooled absorbers next to the tertiary so as to block 
the secondary mirror. We compared the power levels obtained from these 
measurements with those obtained from the receiver's internal calibration 
targets and we obtained Pint = 0.90 ± 0.01. 

Although the reflecting layer, 0.07 pm of aluminum, on the optics is 
only 1/3 skin depth at 183 GHz, calculations following Stratton (1941) 
predicted a reflectivity of 0.994. We placed the loads so as to block the 
tertiary, and thereby determined that the reflection efficiency of the 
tertiary was 0.987 ± 0.005. We thus estimate that the efficiency due to 
reflection at all the optic surfaces, p r , was 0.97 ± 0.01. 
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Since our horn Gaussian beam pattern tapered to -13 db in the E-plane 
and -10 db in the H-plane at the edge of the secondary, we calculate that 
13-14% of the power reflected to the primary was subsequently blocked by 
the secondary (h mo = 0.86-0.87), Thus, the telescope efficiency r )2 = 

0.755 ± 0.01. We also calculated that 6-9% of the power spilled passed the 
secondary (hf S s = 0*925 ± 0.015). 

Figure 2 gives T a (v) calculated using the model of Waters (1976). 

The profile shown is obtained from pressure broadening above 25 km altitude 
(Waters _et_ al^. 1980b). Absorption below 15 km altitude, in that region 
where some variability occurs in the water content of the atmosphere, 
contributes a uniform optical depth of ~ 0.03 across a band exceeding 500 
MHz. Therefore the calculated profile is, for practical purposes, constant 
for all conditions encountered at the KAO operating altitude. An 
interhemispheri c survey found only small variation in the 183 GHz emission 
between latitudes of +60° and -70° (Waters et al . 1980b). The calculated 
average atmospheric opacity across our 100 MHz bandwidth is 0.12. 

The receiver efficiency, rij, and atmospheric transmission calculation 
can be tested against measurements of the atmospheric brightness temperature. 
The antenna temperature observed on "cold” sky can be expressed as 


-T (V )A 
1 - e a 


+ 1 - e 


-T a (Vi)A 


where = 0.755 ± 0.01; T^g is estimated, on the basis of various sensors 
in the telescope cavity, to be 260 ± 10 K; T^jr is estimated from the U. S. 
Standard Atmosphere (1976) to be 250 1 10 K; T a (v s ), the average atmospheric 
opacity in the signal sideband, is 0.12; and T a (vi), the average opacity in 
the image sideband, is calculated to be 0.003. Thus, we calculate that 
TSKY = 78 i 4 K at an airmass of 1.4. We measured the antenna temperature 
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on the sky to be 79 ± 2 K at that airmass. Comparable results had been 
obtained during the 1976 and 1977 flights, indicating that p^ was the same, 
to within 5%, for all the measurements. 

The efficiency factors rj f ss (spillover) and n c (source coupling) were 
calculated assuming a Gaussian field distribution over the primary with no 
scattering at the reflecting surfaces. The half power beamwidth was 
calculated to be 7.9 arcmin in the E-plane and 7.6 arcmin in the H-plane. 

We computed the aperture efficiency, p^ (which can be expressed as 
P r hint hfss times the on-axis gain divided by the on-axis gain for an 
equivalent, uniformly illuminated, unobscured aperture), to be 0.47 ± .02. 

Since there was no point source strong enough to make an observational 
verification practical, we made measurements of the moon. The brightness 
distribution of the moon at 1.64 mm was interpolated from models calculated 
by Keihm (1983) for 1.2 and 3.3 mm, taking into account the orientation of 
the E-vector as projected on the moon. Figure, 3 shows the computed scan 
profile, normalized to the central disk antenna temperature, as well as the 
measured points. The coupling efficiency, p c , for the central position 
was computed to be 0.88, which is about 2% less than that computed for a 
uniformly bright moon, so that we believe that this efficiency may be 
accurate to about - 0.02. The center disk brightness temperature 
interpolated from Keihm's models is 345 ± 10 K. The airmass for these 
observations was 1.4, so that we calculate an antenna temperature difference 
with respect to cold sky of 


AT A T moon ” ^sky 


(345 ± 10 K) 


= *1 ^fss ^c 


-T a (v s )A -T a (Vi)A 


= 195 ± 8 K 
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The measured antenna temperature at the moon's center was 165 ± 3 K. The 
difference between these is about twice our estimated uncertainty and 
represents an additional loss of 15% that we have not accounted for. 

Because the discrepancy is not seen in the sky brightness measurement, we 
suspect our simplified beam pattern calculation. Factors which have been 
overlooked, because they were not convenient to incorporate into our method 
of calculation, include diffraction at the edge of the secondary, and 
obscuration and diffraction by the spider. Based on our actual measured 
antenna temperature on the moon, we therefore derive an overall efficiency 
of 

(hi hf ss n c ) = 0.52 ± 0.05 

and also revise our 1980 aperture efficiency estimate downward by the same 
fraction to n^ = 0.40. 

In 1976 October, the moon was measured to be 153 ± 5 K, corrected for 
atmospheric absorption, at the center of the illuminated portion for which 
a brightness temperature of 315 ± 15 K was estimated (Waters et_ al_. 1980a). 
The lower overall efficiency (nj nf ss P c ) of 0.49 ± 0.05 is entirely 
accounted for by a stronger edge taper (~ 18 db) used to illuminate the 
secondary. We conclude that the efficiencies for extended sources agree 
for the two receivers to within 10%. 

One concern might be the effect of an error in our designed 
illumination of the secondary. Figure 4 shows the computed aperture and 
moon coupling efficiences including the empirical 15% correction. Based 
on the designed edge tapers, we would expect a 15-20% stronger signal from 
an unresolved source in 1980 February than in 1977. Since the difference 
between the designed and actual edge tapers is, at most, one or two dB, we 
can rule this out as a possible cause for the apparent near doubling of the 
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antenna temperature (and certainly would not explain the frequency shift). 
Furthermore, any errors in illumination so gross as to cause this effect 
would also be evident in the moon coupling efficiency. 

c) Results 

Figure 5 shows the observations made in 1977 December (Waters _et_ al . 
1980a) and 1980 February, smoothed to a resolution of 0.5 km s~ 1 . The 
integration times were 110 minutes in 1977 and 19 minutes in 1980. No 
baseline corrections have been applied. The single sideband temperature 
scale, which is twice the antenna temperature scale on the assumption that 
the gains in the two sidebands were equal, has been normalized by dividing 
each data set by the appropriate aperture efficiency. The assumption of 
equal noise power in the sidebands is justified by the narrow sideband 
separation, 2.8 GHz, relative to the observing frequency and the fact that 
we always tuned to the backshort position closest to the diode. 

The difference between the 1977 and 1980 observations cannot be 
attributed to polarization. The 1980 spectrum is an average of data taken 
on two consecutive nights. Differences in the hour angle at the time of 
observation caused the E-vectors to have position angles of 117° and 170°, 
respectively. There were no differences greater than 2 K between the two 
separate observations. Comparison with spectra of other sources showed this 
to be within the noise of the digital spectrometer. If the radiation were 
100% linearly polarized, these two observations could only be reconciled if 
the polarization vector lay in the position angle range of either 5025 to 
56° or 133° to 156°. The receiver E-vector position angle in 1977 December 
varied between 140° and 165°, so that only the 50° - 56° range of polarization 
angles could yield a weaker signal for the 1977 observations. Also, because 
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of variation in the position angle during the course of the observations, 
the intensity of the signal, relative to the 1980 observations, should then 
have varied between at least 50% and 80%, and possibly as much as 0 to 84%, 
depending on the degree of polarization and the position angle of the 
polarization vector. In fact, no such variation was seen during the 1977 
December observations. 

The difference between 1977 and 1980 cannot be attributed to atmospheric 
effects. The vertical dotted line in Figure 5 represents the position of 
the terrestrial H 2 O line in 1977 December. Since the airmass for these 
observations was 1.5, the correction factor at 10 km s’" 1 for losses through 
the atmosphere is 1.19. In 1980 February, the terrestrial line was at 
-39.5 km s' 1 with respect to the local standard of rest; and the airmass 
was 1.6, giving a correction factor of 1.18 at 10 km s'" 1 . The observations 
made in 1976 October and 1977 October (Waters et al . 1980a) agree to within 
the noise, when atmospheric losses are taken into account, with the 1977 
December data. 

A dashed fiducial line has been drawn through the peak of the 1977 
spectrum, at 10 km s _1 , to aid in comparing the spectra. The peak SSB 
antenna temperature in 1977 December was 7 K. In 1980 February, the peak 
was 13 K and had shifted to 11.5 km s'" 1 . Comparing the line wings, we see 
a dramatic increase in the range 15 - 20 km s -1 , but only a little change in 
the range of 0 - 5 km s'" 1 . The upper panel in Figure 5 contains the 
difference between and 1977 and 1980 spectra. The increase at 12 kins'" 1 can 
be clearly seen. The peak intensity of the difference spectrum is comparable 
to that of the 1977 spectrum. The difference profile is narrower than the 
original spectrum, but still has broad wings which extend at least 15 km s'" 1 
on either side of the peak. 
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III. DISCUSSION 


The size of the high velocity source is somewhat dependent on the 
molecule, transition, and technique used to observe it (Table 1). The main 
HVS, at an LSR velocity of ~ 9 km s -1 , appears to be ~ 20 arcsec in 
diameter. If the 183 GHz source were this size, the brightness temperature 
of the additional emission would be ~ 5000 K near the inner part of the 
wings. If this emission arose in only part of the HVS, the brightness 
temperature would be proportionally higher. However, from the shape of the 
profile, we can infer that the emission region is at least large enough to 
include many individual clumps of HV gas. 

In subsection a we present radiative transfer arguments which limit 
the brightness temperature. In subsection b we consider the mechanisms 
which can pump the transition, which also argue against an extremely high 
brightness temperature. We then discuss, in subsection c, possible causes 
of the enhanced emission and, in subsection d,. present the results of 
statistical equilibrium modelling which suggest an explanation for the 
phenomenon. In subsection e we discuss the relation between the 183 GHz 
enhancement and the 22 GHz flare. 

a) Nature of the Emission 

In spite of the change in the spectrum, whiph is usually an indication 
of maser action, the overall appearance suggests that there is little gain 
in the transition. For comparison, consider the 22 GHz H 2 O spectrum and 
the 18 cm OH spectra which consist of numerous spots in the sky emitting in 
very narrow bandwidths (typical 0.5 - 1 km s -1 ). If the rather smooth 183 
GHz spectrum were to be composed of such features, there would need to be 
many more of them. But if the 183 GHz transition should so easily go into 
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a strong inversion, then 183 GHz emission should be quite common in the 
Galaxy. Since it has been seen only in Orion (Waters et_ al_. 1980a), it 
appears that 183 GHz emission is relatively difficult to produce. However, 
the conditions in Orion which favour the inversion of the transition are 
not very selective since the emission occurs without any small-scale 
structure in Doppler velocity and over a wide range of velocities. Since 
velocity coherence does not appear to be necessary to produce the emission, 
we conclude that amplification by stimulated emission does not play an 
important role. 

Since the gain is small, the excitation temperature cannot also be 
small for a strong line to be seen. Thus, it follows that the population 
ratio must be near unity. With the assumptions of small optical depth and 
small population difference, the equation of transfer simplifies to 


3.77 x 10 


-12 Su n u t H 2°l Jh 

n [H 2 ] AV 


where n u is the population per magnetic sublevel of the upper state in the 


the transition, n H 20 is the density of water molecules in all states, g u is 

N 

the statistical weight of the upper level, H 2 is the gas column density, AV 
is the linewidth in km s" 1 , and the other symbols have their usual meanings. 


For illustration, we consider possible emission from the quiescent cloud 
which yields the spike feature. The gas column density is 2 x 10 23 cm*" 2 
(Evans et_ al_. 1975); the linewidth is ~ 5 km s' 1 ; and an upper limit to the 


population in the 3^3 state may be given by 7 exp (-204.5/T)/Q where the 
partition function has been tabulated by Poynter and Pickett (1981). If 
we take 10' 1 * as an upper limit to the water abundance, we find that the 
brightness temperature is < 5 x 10 1 * K. If we assume that the water 


abundance is equal to the CO abundance (Iglesias and Silk 1978; Elitzur 
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1979; Kuiper, Zuckerman, and Rodriguez Kuiper 1981) and estimate the CO 
abundance in Orion to be ~ 2 x 10"® (Wootten eT_ a_l_. 1978), then we find 
that Tg < 1000 K. These crude estimates are consistent with the results 
of statistical equilibrium calculations and show, in a relatively 
straightforward way, that the 183 GHz brightness temperature is quite 
modest compared to those inferred for strong masers (e.g., Reid and Moran 
1981, Elitzur 1982). 

The conditions in the high velocity region are not as well understood. 
The gas column density is significantly less because the CO J = 0 - 1 wings 
are not optically thick (Kuiper, Zuckerman, and Rodriguez Kuiper 1981). 

The abundance of water may be enhanced by shocks (Iglesias and Silk 1978). 
These two effects will tend to cancel each other, so that the conclusion 
that the brightness temperature in the 3^3 - 2 q 2 transition must be modest, 
when compared to prominent astrophysical masers, still holds. 

b) Excitation of the Transition 

We first assess the role of collisions in populating the levels of 
the 3^3 - 220 transition. Figure 6 shows the lower part of the energy 
level diagram for para-H 20 with the allowed radiative transitions which 
connect these levels. The probability of spontaneous decay per second is 
shown next to the lines for the AK = 1 transitions. Collision rates have 
been calculated by Green (1980) to be in the range of 10" 12 - 10 -11 cm 3 
sec -1 for the 3^ and 2 2 q levels over a wide range of temperatures. Thus, 
gas densities on the order of 10 10 cm" 3 would be required in an optically 
thin cloud to excite the rotational levels collisionally . The far-IR H 2 O 
transitions in Orion are, however, optically very thick. Both statistical 
equilibrium and LTE calculations yield estimates in the range of 10 3 to 10 1 * 
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for the extended, quiescent cloud with typical assumptions for temperature, 
density, column density, and water abundance. Consequently, the density 
requirement for collisional excitation is reduced to ~ 10 7 cm"^. In the 
case of clumps of high velocity gas, densities as high as ~ 10 7 cm“ 3 are 
plausible (Kuiper, Zuckerman, and Rodriguez Kuiper 1961; Masson et al . 1983). 
Although the role of collisions may not be entirely negligible, they cannot 
account for the observed high brightness temperature. 

Following Goldreich and Kwan (1974), let us imagine that the 183 GHz 
population ratio is enhanced because of an excess of decays from the 22 q 
level relative to the 3j3 level. (Indeed, the ratios of the appropriate 
probabilities for spontaneous decay, see Figure 6, makes this a plausible 
mechanism.) Therefore, for every 183 GHz photon, there must be at least 
one photon at 2969 GHz. But since the 183 GHz brightness temperature is 
only 10 to 100 times the excitation temperature of far-IR transition, the 
183 GHz emission is enhanced over spontaneous decay by only a factor of 10 
to 100. So, for every 183 GHz photon, there must be approximately 10 3 ” 1 * 3 
THz photons. Now the number of photons emitted per second at 183 GHz per 
unit velocity interval is 

AN/AV = 4ttD 2 (-%■ J • S • (— J/AV 
m \ 4n / \hv/ 

= 10 45 T a (km _1 s ) 

where D is the distance to Orion (460 pc), is 4tt since the gain is small, 

S = 2 kT^/n^A, where A is the physical area of the telescope, and the other 
terms have their usual meanings. Thus, £ 10 4 *® photons per second must be 
emitted by a clump at 3 THz. However, a blackbody at 100 K can only radiate 
lO 45 © 2 photons per second at 3 THz, where 0 is the angular radius of the 
source at the distance of Orion. Thus, only the entire HV source as seen in 
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CO (Table 1) could emit enough far-IR radiation for this mechanism. While 
the emission rate per unit mass increases with decreasing clump size, a lot 
of small clumps in close proximity does not increase the overall loss rate 
from the ensemble because most clumps would radiate into each other. 

The efficiency of collisional pumping by the loss of far-IR line 
photons is enhanced by the presence of cold dust by a factor of the ratio of 
the dust optical depth to the line optical depth (Goldreich and Kwan 1974, 
1979). Gatley et al . (1977) and Harvey, Campbell, and Hoffman (1979) find 
that the absorption efficiency of dust varies as X" 1 from the near-IR to 
the far-IR. Using a ratio of Ay /%2 “ 2 x 10~ 21 mag cm -2 , we find 
t( 102 p) « 10~ 23 N H2* The optical depth in the 22 q In transition is 

1.8 x 10 13 N^(l - n u /n£)/AV, where AV is in km s 1. For an approximate 
evaluation, we assume that the levels are populated according to LTE so that 
x ^20 ~ ln)/ T (102 p) is within a factor of two of 10 3 (%2°/^ H 2)/ AV for the 
temperature range 60 - 200 K. For a predicted abundance of 10~ 6 for water, 
we the escape of photons from the cloud is by several orders of magnitude 
the dominant loss mechanism of 220 “ 111 photons. 

Although the situation is more difficult to analyze for a collisional 
pump which overpopulates the 3^3 level, because of the large number of 
levels involved, the same order of magnitude conclusion applies. Also, as 
we discuss below, the variability argues convincingly against a collisional 
pump. 

Far-infrared radiation, however, does appear to be able to provide the 
necessary excitation. We can use the Einstein relations and the Planck 
distribution to write the ratio of the spontaneous decay rate to the rate 
of radiative excitation as 


- 18 - 



0A 


ill 


8 


(1 - g)B lu WB v (T 1R ) 


(1 ' 3)W 
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where 0 is the probability that a photon will pass through the cloud without 
being reabsorbed, W is a dilution faction which applies to the mean intensity 
of the exciting radiation, incorporating both optical depth and filling 
factor considerations, and the other symbols have their usual meanings. The 
frequencies of the principal transitions by which the 220 ant * ^13 levels 
decay are 3 and 2.2 THz, respectively; so that for temperatures in the range 
of 100 - 150 K, the quantity in square brackets is on the order of unity. 

If the escape probability in the far-IR transitions does not approach unity 
(see above), then the competition between spontaneous decay and radiative 
excitation depends principally on the ratio 0/W. If we assume that the 
source of IR radiation and the emitting gas are generally in close proximity, 
then this ratio becomes the ratio of the far-IR optical depth in the dust 
to the optical depth in the far-IR transitions of H 2 O. As noted before, 
x(100 p) = 10 23 ^ 2 . The optical depths in the 22 q - ljj and 3^ - 2 q 2 
transitions are both 1.8 x 10 13 (1 - n u /n^)/AV. Since the second term 

in the brackets becomes significant only after the upper levels have been 
populated, we may ignore it and note that only one water molecule for every 
10 10 /AV gas molecules need be in the lower rotational states for the 3^3 
and 220 levels to become populated by radiative excitation. If, as we 
previously estimated (Waters et_ al_. 1980a) the abundance of water is ~ 10~ 6 , 
there seems to be no difficulty, for any reasonable velocity dispersion in 
a clump, in having enough opacity in the far-IR rotational lines of water 
for radiative excitation to be effective. 

If, as we have estimated (Kuiper, Zuckerman, and Rodriguez Kuiper 
1981), the abundance of water is higher in the HV gas, then FIR radiative 
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excitation will be even more effective for the HV gas. 

We can easily verify that a far-IR pump can produce enough emission by 
noting that there are sufficient pump photons to produce the observed number 
of 183 GHz photons. On the assumption that the Doppler widths of the 183 
GHz and the pump transition are the same and using the Rayleigh- Jeans 
approximation, one obtains for a saturated maser (e.g., Elitzur 1982) 

T b (183 GHz) < 140 T b (2164 GHz)(ft p /ft 1B 3 ). 
where ftp is the solid angle subtended by the far-IR source at the 183 GHz 
emitting region, and fti83 is the solid angle into which the 183 GHz 
emission is beamed. Since the dust and gas are co-extensive, and the 183 
GHz transition radiates principally by spontaneous emission, these angles 
are both 4 n. We consider first the extended cloud (i.e., "ridge”) for 
which the parameters are better known. Averaged over the central arcminute, 
the emissivity of the dust at 200 - 300 p is ~ 0.1, and the temperature is 
~ 70 K (e.g., Werner et al . 1976) so that T B (-2164 GHz) ~ 5 - 7 K and T B 
(183 GHz) < 1000 K. For a low gain, unsaturated maser, the brightness 
temperature would be significantly less, which is consistent with the 
observations (Waters et al . 1980a). Thus, we conclude that the extended, 
quiescent molecular cloud will emit detectable 183 GHz radiation and will 
have large optical depths in the far-IR water transitions, unless the 
abundance of water is an order of magnitude or more lower than predicted 
by chemical models. Within dense clumps of HV gas, the extinction can 
possibly be higher than the average over one arcminute. Thus, higher 183 
GHz brightness temperatures may be expected for HV clumps. However, it is 
unlikely to be more than one or two orders of magnitude brighter. 
Consequently, the emission region would not be smaller than a few arcseconds. 

We note in passing, the near-IR excitation by hot dust is not likely 
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to be an important mechanism. The model proposed by Goldreich and Kwan 
(1974) is quite general and can be applied with little modification to any 
transition in the ground vibrational state. It involves excitation to the 
first vibration state by means of 6 pm radiation. The mechanism requires 
high density (~ 10 9 cm -3 ) and gas kinetic temperature (~ 1000 K) , and a dust 
temperature several hundred degrees above the gas temperature. However, the 
high velocity gas is much more rarified (£ 10 7 cm~ 3 ) and cooler (~ 100 K) . 

c) Origin of the Emission 

We have previously argued that the high-velocity wings originate in 
clumps of gas from the prestellar molecular cloud which have been accelerated 
to high velocities by the wind of one or more newly formed stars (Kuiper, 
Zuckerman, and Rodriguez Kuiper 1981). Since the dispersal of the original 
molecular cloud is yet in an early stage, the flow is channelled in a mostly 
random and turbulent way through the interstices which have formed in the 
cloud. Thus, the shape of the high-velocity wings does not come primarily 
from an accelerating radially outward flow, but from a high randomized flow 
which has a small net outward component. The significance of this model in 
the current context is that, at angular resolutions insufficient to resolve 
the individual clumps, subdivisions within the region of turbulent flow 
will have roughly the same shape wings and that variations in the mean rest 
velocity among the various subdivisions manifest themselves mainly as 
changes in the velocity of the line center. Random changes of position with 
velocity are seen in high resolution maps in SO (Plambeck et al . 1982), v = 

0 SiO (Wright al_. 1983) and NH 3 (Pauls _et_ al_. 1982). We note that the 
extent of these sources at any velocity greatly exceeds any shift between 
centroids at different velocities (i.e., random motions dominate over 
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systematic ones). The phenomenon also manifests itself in the different 
sizes of the HV sources as seen in various transitions (Table 1). 

Even for a region as small as 10" across, which corresponds to 7 x 10 33 
km, it is hard to imagine any process, other than radiation, which can cause 
a change in three years, since a velocity exceeding 2 x lO 4 km s -3 would be 
required. Thus, one posible explanation for the increased 183 GHz emission 
is that a new region in the Orion complex has become illuminated by 
additional far-IR radiation. 

We can imagine several ways in which the far-IR radiation can be 
increased. These include i) a variable star deeply embedded in the cloud, 
ii) a flare associated with such a star, and iii) a break in the obscuring 
material surrounding such a star. To evaluate these possibilities, it is 
useful to consider the diffusion time for near-IR radiation, which is 
responsible for heating the dust. Using Goldreich and Kwan's (1974) 
estimate for the photon mean free path, the diffusion time in years is 
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We see that at the densities associated with clumps (~ 10 7 cm 3 ), the dust 
could be heated over the entire HV source in 1.6 y. This, of course, is a 
gross overestimate since we believe that the clumps sparsely fill a less 
dense region (Kuiper, Zuckerman, and Rodriguez Kuiper, 1981). 

We consider first the possibility of a change in a dense circumstellar 
obscuring layer being responsible for the change. Yorke and Kriigel (1977) 
have given a model for a dense protostellar cocoon with an inner radius of 
10 33 cm where the gas density is 10 3 ^ cm" 3 , The density varies roughly as 
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r~ 2 . Using the above formula, we see that the radiation takes ~ 10~ 3 y to 
diffuse out through 10* 3 cm, and thereafter doubles its distance from the 
center every 10~ 3 y. Actually, this is an overestimate because as the 
radiation moves outward, the bulk of the energy shifts to longer wavelengths, 
so that the photon mean free path gets longer. Anyway, in about four days, 
the radiation has diffused out 10 17 cm, where the density is ~ 10 6 cm -3 . 
Except in the unlikely case that the density of a cocoon is several orders 
of magnitude higher, the circumstellar dust cannot delay the IR radiation 
long enough to make a physical change in the cocoon a significant event in 
the heating of the outer layers. 

Consequently, any change in the IR radiation must be related to changes 
at the energy source(s). However, the increased 183 GHz emission cannot be 
due merely to an increase in the far-IR radiation which pumps the transition. 
In that case, there would have to be a proportional increase in the far-IR 
luminosity. Gehrz _et_ al^. (1983) have monitored the Orion KL region at 
wavelengths between 2 and 20 ym and have not found any variations larger 
than 30%. Now it is clear from the change in the shape of the profile, as 
well as the probable flare nature of the causal event, that only part of the 
HV region is involved in the change. Thus, a 183 GHz increase due solely to 
increased far-IR pumping, confined to a part of the HV source, would require 
the far-IR increase to be locally even larger, by a factor which is the 
inverse of the ratio of the area of the affected region to the entire HV 
source. Since the amount of dust, and hence the far-IR opacity, can be 
considered constant, a several-fold far-IR flux increase could only be due 
to an increase of the dust temperature to several hundred Kelvin, which 
would certainly be evident near 20 y. Also we would expect the dust to 
cool promptly on a time scale given by Equation 7, which provides the 
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physical basis for the mechanism identified by our modelling calculations. 


d) Statistical Equilibrium Modelling 

We have carried out statistical equilibrium calculations with photon trapping 
to determine the conditions likely to cause a major increase in 183 GHz 
emission. These calculations were carried out. using the large velocity 
gradient (LVG) hypothesis so that radiative coupling between distant parts 
of a cloud can be ignored. In our model, the decoupling arises from the 
large relative separation of the clumps of HV gas as well as differential 
velocities between the clumps. It is not meaningful in this case to use an 
explicit velocity gradient formulation, n H 2 (dV/d£) _1 . Instead, we consider 
explicitly the column density n H 2 dJ l, and the velocity dispersion, AV, in 
the clumps. Sobolev (1957) and White (1977) have suggested that this is a 
reasonable approximation for clouds without velocity gradients. We 
simplified our calculations by neglecting absorption by dust as being 
negligible relative to line absorption (see above), and saturation (i.e., 
the escape probability was forced to be less than or equal to one) because 
high gain cases are not relevant to the present observations. We used the 
collision cross sections computed by Green (1980) for the ten lowest 
rotational levels of para-H20. Because the ~ ^13 transition is one of 
the major pathways feeding the upper level of the 183 GHz transition, we 
extended the published cross sections with unpublished estimates by Green 
(1979) for the rates connecting the next three levels. (The levels 431 and 
5j5 are not shown in Figure 6.) Following White (1977), we expect that the 
numerical results are correct to within about half an order of magnitude. 

The radiation intensity was taken to be blackbody shortward of the 
wavelength, X c , at which the optical depth of the dust is unity. Longward, 
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the intensity was taken as blackbody multiplied by an emissivity which was 
inversely proportional to the wavelength (recall subsection b). It should 
be noted that the quantitative (but not the qualitative) results are 
sensitive to the wavelength dependence of the dust emissivity. Using the 
previously discussed ratio of IR optical depth to total gas column density, 
we computed cases for %2 = 5 x 10 21 , 10 22 , 2.5 x 10 22 , 5 x 10 22 and 10 23 
cm 2 , corresponding to X Q = 5, 10, 25, 50, and 100 pm, respectively. For 
each column density, we computed three cases: gas and dust at 80 K, gas at 

80 K and dust at 200 K, and gas at 200 K with dust at 80 K. (A temperature 
of ~ 200 K is suggested by the excitation temperature of ammonia in the 
"hot core" source — Morris, Palmer, and Zuckerman 1980; Ziurys e_t_ ad_. 1981.) 
In each instance, we carried out the calculations for densities from 10 4 to 
10 9 cm"" 3 by steps of factors of ten. 

We found that, in all cases, a density of 10 9 cm -3 was sufficient to 
thermalize the 3^3 - 2£0 level. With both the gas and dust at 80 K., the 
strength of the 183 GHz inversion decreased with increasing column density 
and, in all cases, was most prominent for a density of 10 6 cm -3 . Brightness 
temperatures significantly in excess of 100 K required that ^i^O/ AV > 10 17 
cm -2 km" 1 s or a relative abundance of water on the order of 10" 5 or higher. 
Increasing the temperature of the dust to 200 K largely suppressed the 
inversion. In this case, the brightest temperatures were ~ 1000 K and 
required ~ 10 18 cm 2 km 1 s, corresponding to a relative water 

abundance of ~ 10 _l *. On the other hand, increasing the temperature of the 
gas to 200 K (with the dust at 80 K) generally increased the brightness 
temperatures. Indeed, if the water abundance is high enough (£, few x 10 5 ), 
it appears that very strong maser emission can occur. Increases as high as 
a factor of 50 were obtained for more modest relative abundances (< IQ- 5 ). 
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For example, at a total gas column density of 2.5 x 10 22 cm" 2 , a relative 
water abundance of 2 x 10"®, a local line width of 1 km s" * , and a density 
of 10 7 (i.e., a source size of 2.5 x 10^® cm or ~ 0.4 arcsec) the brightness 
temperature is 80 K for a gas temperature of 80 K, and 800 K for a gas 
temperature of 200 K. At a density of 10® (i.e., a source size of ~ 4 
arcsec) and a relative abundance of 4 x 10"®, these numbers are 100 K and 
1000 K, respectively. These numbers should be taken only to illustrate the 
qualitative result that an increase in the gas temperature can strongly 
brighten the 183 GHz emission. 

We conclude, therefore, that the observed increase in the 183 GHz 
emission was most likely due to an increase in the temperature of the gas. 

e) Relation to the 22 GHz Flare 

In 1979 October, a 22 GHz 1^0 maser feature at 8 km s 1 flared to reach a 
peak of > 10® Jy (Matveyenko et^ ad. 1980; Abraham et al. 1981). The 
outburst decayed more or less steadily until the middle of 1980 February, at 
which time it was still more than an order of magnitude stronger than any 
previously observed water maser feature in Orion. Late in February, the 
feature flared again, twice as strong as before, and then decayed to a level 
of ~ 10® Jy where it has remained since, varying by ~ 20% (Garay, Moran, and 
Haschick 1983), The general time signature of such events has been shown to 
be consistent with an impulsive injection of energy (Burke, Giuffrida, and 
Haschick 1978). Evidence for impulsive release of mechanical energy 
associated with pre-main-sequence stars is deduced by Mundt and Hartmann 
(1983) from the acceleration of HH1 and HH2, and from the occurrence of 
highly collimated jets (Graham and Elias 1983; Mundt and Fried 1983). 

The 22 GHz flare occurred at a velocity which is in the wing of the 
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183 GHz enhancement, not the velocity at which the 183 GHz line is the 
strongest. Also, we should bear in mind that the conditions which are 
believed to be present in 22 GHz regions would thermal ize the 183 GHz 
transition. The 22 GHz flare region could, however, come from a single 
dense clump in the 183 GHz emission region. This region probably encompasses 
many other 22 GHz masers (Genzel et al . 1981) which have not changed in a 
correlated manner (Cohen and Zisk, 1980). Consequently, any IR radiation 
associated with the 22 GHz flare would have to be at wavelengths short 
enough to be effectively shielded by dust, say, a few microns or less. 

A possible relation between the events at 22 GHz and 183 GHz may be 
less direct than this. For example, mechanical energy associated with a 
flare event could have caused shock heating or possibly a favourable 
alignment of material which, together with IR radiation, could be responsible 
for the 22 GHz maser (Abraham et^ aj^. 1981). The 183 GHz enhancement would 
then be a secondary manifestation of the same flare event in the form of a 
temporary heating of the dust. In a model proposed by Deguchi (1981), a 
burst of infrared luminosity occurs at the center of an optically thick, 
cold, dusty cloud. (In principle, shock heating would serve the model 
equally well.) Once the energy source is turned off, the dust cools more 
rapidly than the gas because the timescales for cooling the gas by radiation 
or collisions with the dust are much longer than the time for radiant 
energy from the dust to diffuse out of the cloud. A hot-gas, cold-dust 
pump, which also selectively removes 523 ' ^14 photons via an absorption 
band at 43 pm, then produces the 22 GHz maser. The 22 GHz emission decays 
in a time set by the rate at which the gas molecules collide with the dust 
which, depending on the density, can be weeks to months (e.g., Goldreich 
and Kwan 1974, Deguchi 1981). 
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The region heated can be considerably larger than the dense (~ 10 9 - 
10 10 cm -2 ) region giving rise to 22 GHz emission. Even at a density of 
10 7 cm -3 , a mean free path at 100 p is only 10* 6 cm (~ 1.5 arcsec), so that 
the photons are considerably more likely to be absorbed than to escape. 

Thus, one would expect, as a natural consequence of Deguchi's model, an 
enhancement of the 183 GHz emission over a larger region. 

There is little difficulty in obtaining enough energy to do the heating, 
To heat the entire plateau source, estimated to be ~ 3 (Zuckerman, 

Kuiper, and Rodriguez Kuiper 1976), from 80K to, say, 200 K would require 
7.5 x 10 93 ergs. Since we require only a small part of the plateau region 
to be heated, let us put the energy requirement to heat the gas at ~ 10 93 
ergs which is only 10" 9 of all the mechanical energy present in the plateau 
(e.g., Bally and Lada 1983). If an impulsive event of this kind occurred 
every twenty years (22 GHz emission was first observed in 1968 by Cheung et 
al . 1969), then the average power required is ~ 5 which is ~ 10 9 of the 
bolometric luminosity of the Orion BN-KL complex, and ~ 2 x 10" 3 of the. 
mechanical luminosity in the high velocity flow. 

At the densities most favourable to strong 183 GHz emission (10 6 ~ 7 
cm" 3 ), the time scale for the heated gas to lose its energy to the cold 
dust is 20 - 200 y (Goldreich and Kwan 1974). If outburst of this kind 

occur as often as every few decades, one might possibly still see evidence 

N 

of previous outbursts. VLA observations of H^ shows most of the hot gas 
to be concentrated in four condensations within a ~ 10” diameter region near 
the center of the Orion IR complex (Pauls et al . 1983). These condensations 
are not responsible for the 183 GHz enhancement or the 22 GHz outburst 
because of discrepancies between the LRS velocities. Also, the Hot Core was 
recognized in observations as early as 1977 (Morris, Palmer, and Zuckerman 
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1980). Townes e:t_ al . (1983), who observed the a(4,3) - s(3,3) transition of 
NH 3 at 125 pm from the hot core, point out however that the gas is at least 
several tens of degrees hotter than the associated dust. Since only one of 
these condensations is associated with a known infrared source (IRc 4) and 
it is questionable whether that one has an internal energy source (Werner 
et al. 1983), the gas will eventually lose its energy by collisions with the 
dust unless it is being continuously resupplied with energy. Thus, it is 
conceivable that these hot spots are earlier manifestations of the kind of 
event which we infer to have occurred in 1979. 

On a speculative note, if the flare event seen at 22 and 183 GHz is a 
by-product of a kinetic event which released considerably more energy than 
is manifested in the heating of gas, then it is possible that such sporadic 
outbursts provide the energy and momentum for the high velocity gas. 

IV. CONCLUSIONS 

The peak flux of the 183 GHz water line emission from Orion increased 
by 50% between 1977 December and 1980 February. An increase was also seen 
in the wide wings of the line. The profile shifted approximately 1.5 km s~* 
to the red. We considered instrumental effects, absorption in the earth's 
atmosphere, and polarization effects as possible causes for the observed 
increase and have been unable to identify any sources of systematic error. 
Simultaneous, undetected errors in both the frequency and intensity 
calibrations seem to us very unlikely. 

Because of the change in the shape of the profile, we concluded that 
the additional emission came from a subset of the entire high velocity 
("plateau”) source. The wide range of velocities over which the increase 
was seen suggests to us that it is not due to a high gain maser, but the 
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Inferred brightness temperatures are high enough to require low-gain maser 
action. Collisions do not appear to be capable of providing sufficient 
excitation in HV clumps because they cannot rid themselves of far-lR photons 
fast enough to support the pump. On the other hand, there is enough far-IR 
emission from the dust to excite the rotational ladder of water and produce 
the observed brightness temperatures. We carried out statistical equilibrium 
calculations which indicate that the likely cause of a large increase of the 
183 GHz emission is impulsive heating of the gas and dust, but with the dust 
having quickly returned to approximately its original temperature. 

The occurrence of a major flare in the 22 GHz emission during the 
period of the 183 GHz brightening suggests a common cause. Because of the 
high velocity of propagation which would be required for mechanical (shock) 
heating of the 183 GHz region, a radiant heating seems likely, although the 
ultimate cause of the radiation may still be a a mechanical event. The 
energy required to cause the heating is only lO -4 of all the mechanical 
energy in the HV source. Even if such an event occurred only every few 
decades, the average power required is ~ 10 _t * of the bolometric luminosity 
and ~ 2 x 10~ 3 of the mechanical luminosity in the Orion BN-KL complex. 

Separate evidence for hot gas in the presence of cool dust is found in 
the hot core source, studied principally in ammonia. The components of this 
source may be remnants of earlier flare events. The timescale for cooling 
of the gas responsible for 183 GHz emission by collisions with the dust is 
expected to be in the range of 20 to 200 years. We also point out that if 
the heating of the gas, which is the proposed cause of the 183 GHz flare, is 
a by-product of a causative kinetic outburst, then the energy and momentum 
in the outburst might be large enough so that a succession of outbursts 
could be the driving force for the high velocity gas. 
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Table 1 - Sizes of HV Source in Various Transitions 


Molecule 

Transition 

Size 

Xlsr . 

Reference 



arcsec 

km s _1 


H 2 O 

6 1,6 ' 5 2,3 

30" x 45" 

~ 10 

Genzel et al . 1981 

CO 

1-0 

30 - 40 

9.5 

Solomon et al_. 1981 


1 - 0 

~ 20" 

9 

Masson _et_ al . 1983 



~ 11" 

5 

Masson et al. 1983 


2 - 1 

~ 25" 

9 

Knapp et_ al . 1981 


3-2 

30 - 40" 

9 

Erickson et a\_. 1982 

HCO + 

1 - 0 

< 30" (1) 

~ 9 . 

Olofsson et al . 1983 

SO 2 

8 1,7 ~ 8 0,8 

~ 25" 

8.3 

Schloerb et_ a 1 . 1983 


13 1 , 13 - 12 0,12 

~ 25" 

9 

Knapp et_ al . 1981 

SO 

2i - 1 X 

11" x 9" 


Plambeck et al. 1982 

nh 3 

(3, 3) 

10" x 7" 

~ 5 

Bastein et al_. 1981 





Pauls et_ al . 1982 

SiO 

1 - 0, v = 0 

< 10" 

~ 5 

Genzel et_ jrl.. 1980 


2 - 1 , v = 0 

15" x 30" 

5-9 

Wright et al. 1983 

(1) With 

a medium velocity 

feature extending 

; 100" to 

the north. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


- Schematic diagram of receiver and telescope system for the 
1980 August flights. As noted in the text, there was a minor 
difference in the method of calibration during the 1980 February 
flights. 

- Optical depth through the atmosphere towards the zenith as a 
function of frequency relative to 183310.09 MHz. The horizontal 
line representing a range of Doppler velocities from -50 to +50 
km s" 1 is provided to facilitate comparison with the Orion 
spectra in Fig. 3. 

- Measurements of a scan across the moon's equator are shown 
along with the results of convolving our calculated beam profile 
with a model of the lunar brightness temperature at 1.6 mm, as 
discussed in the text. 

- Aperture efficiency and moon coupling efficiency as a function 
of the illumination at the edge of the secondary reflector 
(edge taper). These efficiencies include the empirical 15% 
correction discussed in the text. 

- Observations of the 3^3 - 220 transition of water in the direc- 
tion of the Orion Nebula in 1977 December and 1980 February. 

It is assumed that the SSB antenna temperature is twice the 
observed (DSB) antenna temperature. 

- The lower part of the para-^O rotational energy ladder, show- 
ing the allowed transitions which connect the levels, and the 
Einstein A coefficients (probability of spontaneous decay per 
second) for the AK = 1 transitions. 
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